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An Isotopic Fingerprint of Electron-Phonon Coupling in High-Tc Cuprates 
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Angle-resolved photoemission spectroscopy with low-energy tunable photons along the nodal di- 
rection of oxygen isotope substituted Bi2Sr2CaCu2 0g+i reveals a distinct oxygen isotope shift near 
the electron-boson coupling "kink" in the electronic dispersion. The magnitude (a few meV) and 
direction of the kink shift are as expected due to the measured isotopic shift of phonon frequency, 
which are also in agreement with theoretical expectations. This demonstrates the participation of 
the phonons as dominant players, as well as pinpointing the most relevant of the phonon branches. 

PACS numbers: 74.70.Pq, 74.25. Jb, 79.60.-1 



Effects of electron-boson interactions on electronic self- 
energies show up as sudden changes in the electron dis- 
persion, or "kinks" , via angle-resolved photoemission 
spectroscopy (ARPES) While the kinks are now 

indisputable, their origin as arising from electronic cou- 
pling to phonons [IjHiQ, magnetic excitations 
or bothflO], remains unclear. The best way to iden- 
tify the origin of an interaction is to slightly modify the 
interaction in a controlled manner. For phonons, this 
can be done via an isotopic exchange, which varies par- 
ticle masses and hence the vibrational energies. Here 
we substitute the oxygen isotopes (^^O ^^O) lead- 
ing to a softening of phonon energies of a few meV 
[=Q X (1— vlG/lS)] where are the bare phonon fre- 
quencies of 40-70 meV [ill [l3|- Recent scanning tun- 
nelling microscopy (STM) experiments showed an iso- 
topic shift of a second derivative feature of 3.7 meV 
[l3|, though whether the observed shift is evidence for 
electron-phonon coupling has been disputed by a num- 
ber of g roups in terms of the inelastic tunnelling barrier 



lj,[15|,[16|. Even earlier than that, Gweon et al. used 
ARPES to study isotope effects on the electronic struc- 
ture _l3l. They found that the primary effect was at 
high energy and was unusually strong (up to 30 meV), 
though these results were not reproduced by more mod- 
ern experiments £1, Therefore, a direct and clear 
"fingerprint" of strong electron-phonon interactions in 
the high-Tc superconductors has been missing up to now. 

The momentum-selectivity of ARPES brings an addi- 
tional tool to bear on the study of the isotope effect which 
is not available from tunneling. This information is very 
helpful in separating out the energy and momenta of the 
bosonic modes which couple most strongly to the elec- 



tronic degrees of freedom. Here we focus on the "nodal" 
spectra, i.e. those along the (0, 0)-(7r, tt) direction in 
the Brillouin zone where both the superconducting gap 
and pseudogap are minimal or zero [H, This is the 
portion of the Brillouin zone where the kinks or self en- 
ergy effects have been most strongly studied both exper- 
imentally 0, S HH and theoretically S H^. As de- 
scribed very recently by Giustino et al. [2J] and Held et 
al. [i^, there is still great controversy about the nature 
of the coupling in this direction; these two theoretical ar- 
guments imply that the phonon coupling could only be 
a minor player in the overall electron-boson coupling of 
these states because the calculations indicate a signifi- 
cantly weaker kink than what is found in experiment. 

Analysis and interpretation of the nodal data is also 
more straightforward than the data throughout the Bril- 
louin zone, as the energies of any bosonic modes in the 
spectra are expected to be shifted by the energies of 
the gaps. Assuming any isotope shifts go roughly as 
[n X (1— VTG/IS)], these are at most a few meV, requir- 
ing an unusually high precision ARPES experiment. To 
verify such a tiny effect in the most convincing way, we 
utilized bulk-sensitive low photon energy ARPES (LE- 
ARPES), which has greatly improved spectral resolution 
compared to "conventional" ARPES [26|. In this Letter, 
we report an isotopic fingerprint of electron-phonon cou- 
pling as a clear isotope shift of the kink energy taking ad- 
vantages of LE- ARPES. Present results provide straight- 
forwardly that the origin of the nodal kink is due to the 
electron-phonon interactions. 

High-quality optimally doped Bi2Sr2CaCu2 08+5 
(Bi2212) single crystals were prepared by the travelling- 
solvent floating-zone technique [131 ■ Oxygen isotope 
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FIG. 1: (color online), (a) and (c) Fermi surface maps taken with 7.5 eV and 7.0 eV photons, respectively, (b), (d), (f) 
Energy-momentum distribution maps taken with 7.5 eV, 7.0 eV and 24.0 eV along the nodal direction, indicated by the red 
line in (a) and (c). MDCs at the Fermi level {Ef) axe also shown by the white line, where the BB and AB denote bonding band 
and antibonding band, respectively. The Fermi momentum is set to the momentum of the AB crossing Ef- (e) Comparison of 
MDCs at Ef between 7.0 eV (blue) and 24.0 eV (red). Note that ARPES data using a conventional photon-energy of 24.0 eV 
were taken at BL5-4 of the Stanford Synchrotron Radiation Laboratory (SSRL). 



substitution was performed by annealing procedures 
[2^ . yielding a slight decrease of from 92.1 K to 
91.1 K. A softening of the oxygen vibration modes 
and the high oxygen isotope substitution rate (more 
than 80 %) was confirmed by Raman spectroscopy [2^. 
Further, LE-ARPES gives us unparalleled precision in 
determining the Fermi surface areas directly from the 
Fermi surface maps, and we have confirmed that these 
are the optimal and the same for the two isotopes to the 
level below 2 % (a doping level uncertainty Z\x=.003). 

Present data were collected at BL-9A of Hiroshima 
Synchrotron Radiation Center using a Scienta R4000 
electron analyzer. We used our newly developed high- 
precision 6-axis sample manipulator to remove the ex- 
trinsic effect due to a sample misalignment. The clean 
and flat surface of the samples was obtained by cleav- 
ing in situ in ultrahigh vacuum better than 4 x 10~^^ 
Torr below 10 K. The total instrumental energy and an- 
gular (momentum) resolutions were better than 5 meV 

and 0.4° (0.005 A~^) at 7.0 eV photons, respectively. 

First, we optimized the excitation energy to avoid 
the complications in extracting accurate band disper- 
sions due to the very small but finite (0.01 A ) bilayer- 
splitting which is now known to exist even along the 
nodal line [Figs. [T] (a) and (b)] js^l- By tuning the pho- 
ton energy to 7.0 eV, we can see the complete isolated 
and individual antibonding band dispersion [Figs. [T] (c) 



and (d)], giving the extremely fine quasiparticle disper- 
sion with the momentum full width (Ak) below 0.005 

A [blue line in Fig.[l](e)]. This is in sharp contrast to 
the cases of the dispersion with bilayer-splitting [Fig. [1] 
(b)], and also to the conventional ARPES spectra broad- 
ened by a limit of resolutions as well as the unresolved 
bilayer-sphtting [Fig. [T] (f)]. These advantages of LE- 
ARPES, including the band selectivity, are mandatory 
to verify the subtle change in the ARPES spectra with 
isotope substitution. 

Figure O (a) compares the ^^O (blue) and ^^0 (red) 
energy-momentum dispersions near the nodal region 
taken with 7.0 eV photons. We found a very small iso- 
tope effect around the kink in the dispersions, indepen- 
dent of the cut position. We emphasize here that each 
of the five positive dispersions is exactly identical with 
each of the five negative dispersions for both ^^O and 
[as seen in Fig. [2] (b)], which rules out the possibil- 
ity of an extrinsic effect due to a sample misalignment. 
An earlier attempt at studying isotopic shifts by conven- 
tional ARPES showed many unusual effects of the oxygen 
substitution, characterized by extremely large (30 meV 
scale) shifts of the high binding energy portions of the 
spectra [l3|- These unusual isotope effects were not ver- 
ified in recent ARPES experiments 18, nor by the 
ultra-high resolution experiments shown here. It should 
be noted that the present results are consistent with the 



3 




150 100 50 

Binding energy (meV) 



FIG. 2: (color online), (a) Energy-momentum dispersions 
near the nodal region both for ^''O (blue) and (red) from 
optimally doped Bi2212. Those measured cuts are labeled 
-5 to +5, displayed in the inset. Gray shaded are roughly 
indicate the kink in the all dispersions, (b) Real-parts of 
the self-energy both for ^^O (blue) and ^^O (red), showing a 
isotope shift of ~ 70 meV peak. Blue and red dashed line 
indicates the kink energy, averaged from cut -5 to cut -f5, 
for ^^O and ^®0, respectively. Positive and negative cuts are 
represented by solid and dashed lines, respectively. 



previous our ARPES [3, [IS] , though wc can not see a 
convincing evidence of such a tiny isotope effect in there. 
This should be attributed to that the experimental reso- 
lutions and/or accuracy were still not enough to observe 
a few meV order effect. 

To visualize the subtle isotope effect more clearly, we 
deduced the real part of the self-energy ReIJ{uj). We ex- 
tracted KeE^Lu) by subtracting a bare band dispersion 
from the experimental one using three different forms of 
bare band in this analysis: linear (presented in the fig- 
ure) [sH, second order polynomial [3^ . and third order 
polynomial. An important thing is that the present re- 
sult is robust in regard to this choice. As seen in Fig. [2] 
(b) , clearly observed is a few meV isotopic shift of the ~ 
70 meV peak of Rei7(a;). Here we estimated the kink en- 
ergy as the energy giving the peak maximum of HelJiu;), 
which was obtained by fitting the top portion of IieS{uj) 
with a Gaussian. To quantify the energy-scale of the kink 



as well as the isotopic kink shift, we have studied multiple 
samples systematically. 

Figure [3] (a) shows the real part of the self-energy 
ReSiui) both for ^^O (blue) and ^^O (red) multiple sam- 
ples. From these Rei7(a;), we obtained kink energy plot- 
ted as a function of five different samples each for both 
^^O and ^^O in Fig.[3](c). We found a clear isotopic soft- 
ening of the kink energy from about 69.0 meV to about 
65.6 meV, or a softening of 3.4 ± 0.5 meV [s^. Addition- 
ally, we used a completely independent analysis method 
using the widths of the ARPES peaks. This analysis 
has the advantage of not having any assumptions about 
a bare band, such that the isotope effect should appear 
more straightforward. Thus we see a 3.2 ± 0.6 meV shift 
in the imaginary part of the self-energy lmE{u}) [Figs. [3] 
(b) and (d)]. By studying ten samples as well as by using 
multiple independent analysis methods, we compensated 
for possible systematic errors which might come into play 
when trying to determine energies to such a great preci- 
sion. Therefore, we can state with confidence that the ~ 
70 meV feature in the nodal electron self-energy is due 
to the coupling of the electrons with phonons. That this 
is the dominant feature in the electron self-energy, as is 
seen from both the real [Fig. [3] (a)] and imaginary [Fig. [3] 
(b)] parts of the spectrum, is clear and significant. 

Then, which phonons are responsible for this cou- 
pling? Neutron scattering experiments 3j| as well as 
first-principles phonon calculations (35l | indicate a few 
phonon modes that are likely to be most relevant for the 
coupling; the in-plane "half-breathing" phonon mode (Jl 
^ 70 meV) and the "buckling/stretching" modes (il ^ 
36 meV). The apical oxygen stretching mode (fi ^ 50 
meV) could also be considered, though calculations in- 
dicate that the number of allowed final states for these 
phonons is negligible 2J|. By coupling the nodal elec- 
trons with momentum, k, to other parts of the Fermi 
surface (k'), the electron self-energy can in principle 
pick up the energy of the superconducting gap /i(k') 
weighted over the Brillouin zone (and hence gap size) by 
the electron-phonon matrix elements |(/(k, k')|, though 
there is still theoretical disagreement about how to do 
this gap-referencing in a d-wave superconductor [IHj ■ On 
the one hand, it is suggested that the kink energy ap- 
pears at the mode energy plus the maximum gap energy, 
or fl+Amax (ill- This scenario would indicate that the 
nodal kink is caused by the 36 meV buckling phonon. Al- 
ternatively, it is suggested that the kink represents sim- 
ply the mode energy, indicating that the half-breathing 
phonon is dominant and couples electrons strongly along 
the node [IS, HE HI] ■ It is here that the present isotope 
effect is uniquely capable to distinguish the two. 

For the buckling mode with ~ 36 meV, the isotope 
shift would be expected to be 2.1 meV, more than two 
error bars outside of the results published here. On the 
other hand, the breathing mode with ~ 69 meV would 
have an isotope shift of 3.9 meV, within the error bar 
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FIG. 3: (color online), (a) Real part of the self-energy ReE{Lo) from five samples both for ^''O (blue) and ^**0 (red) along the 
nodal direction indicated by the red line in Fig.[T](c). All Rei7(tj) are deduced by subtracting a bare band dispersion from the 
experimental one, where w is the energy relative to the Fermi energy, and normalized by the peak maximum, and are also offset 
for clarity, (b) Imaginary part of the self energy Im£'(a;) determined from MDC full widths. An impurity scattering term at 
ci; = is subtracted as an energy independent constant background, (c) and (d) Obtained kink energy as a function of sample 
numbers both for ^''O (blue) and (red) from ReX'(a;) and lmE(ui), respectively. 



of the isotopic shift seen here. This serves as a strong 
indication that the breathing mode is responsible for the 
nodal kink in ARPES data, consistent with expectations 
from recent theory 2J] . However, the moderately strong 
coupling that we observe, (a coupling parameter A ^ 0.6 
[3^) is significantly stronger than that calculated by the- 
ory [13, [2I], indicating that these calculations are missing 
an important ingredient to the electron-phonon coupling 
in the cuprates. Nominally we expect such enhancements 
in the electron-phonon coupling to originate in the strong 
electron correlation effects — a problem which is in- 
teresting in its own right, but which also may have par- 
ticular relevance to the mechanism of superconductivity 
in the cuprates. 

In summary, we reported the isotopic fingerprint of the 
nodal kink probed by high-precision LE- ARPES. Present 
isotope shift of the kink energy provides the first convinc- 
ing and direct evidence that the electron-phonon inter- 
actions are responsible for the origin of the nodal kink. 
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